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ABSTRACT 

The energy spectra of ultra-high energy cosmic rays (CRs) measured with giant extensive air shower 
(EAS) arrays exhibit discrepancies between the flux intensities and/or estimated CR energies exceeding 
experimental errors. The well-known intensity correction factor due to the dispersion of the measured 
quantity in the presence of a rapidly falling energy spectrum is insufficient to explain the divergence. 
Another source of systematic energy determination error is proposed concerning the charged particle 
density measured with the surface arrays, which arises due to simplifications (namely, the superposition 
approximation) in nucleus-nucleus interaction description applied to the shower modeling. Making 
use of the essential correction factors results in congruous CR energy spectra within experimental 
errors. Residual differences in the energy scales of giant arrays can be attributed to the actual overall 
accuracy of the EAS detection technique used. CR acceleration and propagation model simulations 
using the dip and ankle scenarios of the transition from galactic to extragalactic CR components are 
in agreement with the combined energy spectrum observed with EAS arrays. 
Subject headings: cosmic rays - instrumentation: detectors - methods: data analysis 



1. INTRODUCTION 

Ultra-high energy cosmic rays (UHECRs) are presently 
measured using a number of giant extensive air shower 
(EAS) arrays. The observed energy s p ectrum ex- 
hibits the cutoff predicte d by iGreisenl (|1966f ) and 
iZatsepin and Kuzminl (|1966f) (GZK) and 'Ankle' features 
at energies of app roximately 4 x 10 19 eV and 5 x 10 18 eV 
(jFukus hima 2009). However, there are essential discrep- 
ancies between the flux intensities and/or the estimated 
energies of the initial cosmic rays (CRs) generating EASs 
detected by different arrays. These discrepancies exceed 
instrumental errors and it makes it difficult to decide for 
the validity of the results obtained. 

The current paper presents an analysis of the sources 
of these discrepancies. It is shown by i) correcting the 
CR intensity due to instrumental errors and power law 
spectrum and ii) taking into account model uncertainty 
in the estimation of EAS initial nucleus energy that the 
observed UHECR energy spectra appear to be congruent. 
Residual differences in UHECR energy scales of arrays 
can be attributed to the actual overall accuracy of the 
EAS detection technique. 

2. TECHNIQUES USED TO MEASURE THE 
ENERGY SPECTRUM OF UHECRS 

There are two basic techniques for measuring UHECR 
parameters with EAS arrays and, in particular, for esti- 
mating the energy. The first is the measurement of the 
electromagnetic component and/or muons reaching the 
ground; the energy is estimated using a model simulation 
of the particle density at a particular distance from the 
shower axis (e.g., Sgoo)- The secon< i technique is based on 
the measurement of the ionization integral of the longi- 
tudinal EAS profile with fluorescence or Cherenkov light 
detectors. In this case the UHECR energy, Eq, is esti- 
mated as a sum of the ionization integral, £j, and the 
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'missing energy', E m , carried by hadrons, muons, and 
neutrinos. The ionization integral is given by 



where N e (t) is the number of electrons and positrons at 
depth t; e is the critical energy in air; to is the electron 
radiation length in air. The missing energy is compar- 
atively small (E m /Eo < 0.1) at energies above 1 EeV 
(=10 18 eV), so the method used can be considered to be 
nearly model-independent. 

A typical example of an instrument applying the 
first technique i s the Akeno Giant A ir Shower Array 
(AGASA) (e.g., iTakeda et all (|200l ). while the sec- 
ond approach is realized in the High Resolution Fly's 
Eye (HiRes) array, consis ting of fluorescence light detec- 
tors (jAbbasi et al.l (|2008l ) and references therein). Next- 
generation arrays combine both techniques : The Pierre 
Auger Observatory (PAO: iSchussler et al.l (120 09)) and 
the Telescope Array (TA; iBergman et al.l (1200 9)) com- 
prise charged particle detectors and fluorescence tele- 
scopes. 

In the Yakutsk array experiment these two techniques 
are also realized: There are scintillators on the ground 
detecting electrons, positrons, photons, and muons; scin- 
tillators beneath the ground detecting muons; and pho- 
tomultiplier tubes (PMTs) detect ing the air Cherenko v 
light produced by the showers (jEgorova et al.l [2004). 
Data from the scintillators and PMTs are used to esti- 
mate the energy of the UHECR particles initiating EASs. 
It has been shown that the two independent estimates 
diverge significantly, the shift in the lg E between resul- 
tant spectra o f CRs is approxima tely 0.12 in the range 
E n > 10 18 eV (<Ivanov et al.|[2009l ). 

3. A CORRECTION TO THE CR INTENSITY 
MEASURED WITH EAS ARRAYS 
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Fig. 1— The differential UHECR flux (multiplied by E$ and 
correction factors from TablefTt measured wi th EAS arrays . Expe r- 
imental data are from AGASA (rhombuses, Takcda ct al. (2003)), 
HiRes I, II (open and fi l led tr iangles, lATibasi ct al. ( 2008j)), TA 
(crosses, Bergman ct al. (2009)), Haver ah Park (HP, h o rizont al 
bars, lAve et al.1 (120031 1), PAO (squares, ISchussler eTaH ll200li)). 
the Y akutsk array (open, crossed and filled circles. lEgorova et al.l 
12001 )). 



A first step in comparing the observed energy spectra is 
to include a correction to the measured intensity of CRs 
caused by the inst rumental errors and power spectrum. 

It was shown bv lZatsepinl (|1959( ) that there should be a 
difference between the observed intensity of CRs and the 
original intensity in the case of a rapidly falling energy 
spectrum, due to instrumental errors and fluctuations in 
the sho wer parameters mea sured. 

Then IMurzml ()1965D and IKalmvkovl (|1969D calculated 
the measured intensity in the case of a lognormal distri- 
bution of Seoo and the so-called shower size, iV e : 

J(N e ) = J (N e ) exp( N \ 2 N) ), (1) 

where Jo is the actual intensity; <7jv is the RMS deviation 
of lnN e ; k is the integral energy spectrum index; and 
ajy = d In N e /d In Eq . 

In our case, the target values are the parameter E = 
'primary particle energy' estimated after shower detec- 
tion, and the actual energy of the CR, Eq, that initiated 
the EAS. The estimated energy has a distribution around 
the mean value formed by the instrumental errors and 
fluctuations with a RMS deviation, a. The energy fluc- 
tuation is small in comparison with instrumental errors 
and can therefore be neglected. 

If we assume the lognormal distribution of y = In E, 
with an average value equal to In Eq , then the observed 
intensity of CRs is given by the convolution of the pri- 
mary spectrum, Jrjexp(— nz), and the distribution of in- 
strumental errors: 

f°° exv(——) 
J(z) = Jo / exp(-Kz + ny) * K 2 ° 2j dy. (2) 

The re sultant observed-to- initial intensi ty conversion fac- 
tor is (|Murzinlll965tlIvanov et al.ll2009h : 

a 2 K 2 

Rj = Mz)/J(z) = exp( — ). 



TABLE 1 

Instrumental errors, ct, and intensity correction factors, 

Rj, FOR EAS ARRAYS. 



Array 


AGASA 


HiRes 


TA 


Yakutsk 


CT, % 


25 


23 


23 


32 


Rj(E < 4.47 EeV) 


0.86 


0.88 


0.88 


0.77 


Rj(AA7 < E < 56.2 EeV) 


0.90 


0.92 


0.92 


0.85 


Rj(E > 56.2 EeV) 


0.53 


0.59 


0.59 


0.36 



The necessary conditions are a constant index and RMS 
error. As a crude approach, one can use the broken power 
law approximation of the energy spectru m given by the 
HiRes collaboration (Abb asi et al.l |20 0Sf) . and the con- 
stant RMS error averaged in the range Eq > 10 18 eV. 
In the vicinity (with width a) of the break points the 
interpolation of the index can be used to prevent gaps 
in the spectrum. The procedure is inevitably iterative: 
The revised spectrum indexes (2.24, 1.8, and 4.5 below 
the ankle, between the ankle and the cutoff, and above 
the cutoff energies, respectively) rather than observed 
indexes should be used in a correction factor. 

In Fig. [T] the energy spectra presented were observed 
by giant EAS arrays with CR intensity correction factors 
calculated using the HiRes power law approximation of 
the energy spectrum and instrumental errors inherent 
to arrays (Table Q]). In the PAO and HP cases, how- 
ever, t he corrections are already applied in t he original 
works (|Ave et al.ll200l ISchussler et al.l I2009Q so no in- 
tensity corrections are made. 

Concerning the Hi Res experiment, the monocular re- 
construction results (jAbbasi are used here 
to yield spectra with the best statistical power over a 
wide energy range. Energy estimation errors of the two 
HiRes detectors in the monocular mode are derived based 
on the original data (jAbbasi et al.lr2009al) . with the ra- 
tio distribution of energies measured by HR1 and HR2 
independently for the same EAS event. The RMS devi- 
ation of the ratio is found to be 0.33 ± 0.01. An imme- 
diate consequence is the average energy estimation accu- 
racy of the Fly's Eye detectors in the monocular mode, 
SE/E ~ ^0MS' 2 HR i +S 2 HR2 ) = 0.23 ± 0.01. 

Sinc e the preliminary r esults from the TA experi- 
ment ((Bergman et al. 2009) are obtained by fluorescence 
detectors in monocular mode applying the same data- 
handling procedure as in the HiRes case, the same en- 
ergy estimation errors are assigned here to both arrays. 
For other experiments, instrumental errors in energy es- 
timation are taken from the original papers. 

4. ENERGY SCALE DIFFERENCE DUE TO THE 
EAS MODELING UNCERTAINTY 

The energy of the primary CR particle initiating the 
EAS is estimated using model relations between Eq and 
measured shower parameters, such as the ionization inte- 
gral, Ei, the particle density at 600 m from the axis, ^ooj 
and the number of electrons at observational level, A^ e . 
These relations are more or less dependent on the high- 
energy hadron-nucleus and nucleus-nucleus interaction 
models used. Thus, in addition to including instrumen- 
tal errors, we must introduce a systematic 'EAS modeling 
uncertainty' into the energy estimation procedure. 

Existing inconsistencies in the interpretation of the ex- 
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Fig. 2. — Probability distribution of the wounded (interacting) 
nucleons fraction, n w /A, in eikonal approximation. Calculation 
results for collisions of nitrogen and iron nuclei in air are shown. 
Average values are indicated by the vertical bars on the bottom 
scale. 

perimental data on the primary energy and mass com- 
position of UHECRs indicate the presence of this uncer- 
tainty. A number of attempts have been made to esti- 
mate the value of the effect and elucidate th e source(s); a 
typica l example is an analysis conducted bv lKnapp et al.l 
(2003 ) of EAS si mulat ions at energies above 10 EeV. 

IKnapp et al.l (|2003f ) found a clear trend of conver- 
gence between different hadronic interaction models: 
The ionization integral differs by only 2% in QGSJET 
01, SIBYLL 2.1, and DPMJET 2.5, but the lateral dis- 
tribution of EAS particles far from the core (r > 600 m) 
still differs by about 15% (photon and electron densities) 
and 30% (muons) in the models. This means that there 
is an uncertainty in the energy determination at the sur- 
face water-Cherenkov detectors (HP, PAO) of about 20%, 
and at the scintillation detectors (AGASA, Yakutsk) of 
~ 15% for vertical showers and up to 30% for inclined 
showers where muons dominate the charged particle den- 
sity at zenith angles 9 > 60°. 

This estimate is undoubtedly a lower limit of the ac- 
tual uncertainty, because it does not include the contri- 
bution from low-energy hadronic interacti ons and com- 
putati onal and coding errors. For instance, IKnapp et al.l 
( 2003) have compared the performance of CORSIKA and 
AIRES codes for the same interaction model (QGSJET). 
The resultant mean positions of EAS maximum in atmo- 
sphere differ by about 25 g cm -2 , the electron numbers 
at the maximum differ by 6%, while the photon and elec- 
tron densities at 10 3 m from the core differ by ~ 20%. 

All these sources of energy estimation errors are sym- 
metric in the sense that they can equiprobably increase 
or decrease the derived primary energy due to an inade- 
quate interaction model or coding errors. 

However, there are other examples of the energy di- 
vergence sources that result in the systematic errors of 
E. The first source is the variation in the mass compo- 
sition of UHECRs. In this case we can use CORSIKA 
simulation results for the primary proton and iron nu- 
cleus (leaving aside more exotic primaries such as pho- 
tons, neutrinos, etc.), QGSJET and SIBYLL models, 
applied to AGASA, Yakutsk array scintillation detec- 



FlG. 3. — Relative difference in charged-particle density, S600 
(circles), and ionization integral, Ei (squares), in fragmentation 
and superposition models. 

tors (jNagano et alj|2000t IDedenko et al.ll2009D and PAO 
fluorescence-detectors (jBarbosa et al.l 120041 ) In the case 
of iron primaries the energy estimate should be decreased 
by ~ 10% (AGASA), ~ 15% (Yakutsk) and ~ 3% 
(PAO) in comparison with proton-initiated showers. A 
corresponding co rrection in the HP case is applied by 
lAve et all (l200l . 

Another example is the use of a common superposition 
model when treating nucleus-nucleus collisions in EAS 
simulation codes. While in some models (e.g. QGSJET, 
DPMJET) the variants of the nuclei fragmentation are 
implemented, in others (SIBYLL, HDPM) only the su- 
perposition approximation is used (jHeck et al. 1998|). In 
contrast, real nucleus-nucleus interactions are mostly pe- 
ripheral, and the probability that the 'wounded' (inter- 
acting) nucleon number, n Wl is equal to the projectile 
mass number A (as in superposition approximation) , is 
small. Fragmenting primary nucleus induces a shower of 
secondaries slowed down in comparison with EAS in the 

superposition model. This should result in S[ Q0 > Sq 00 
and then a reduced primary energy estimate in the frag- 
mentation model in comparison with the superposition 
model. 

In the next section, EAS observables are estimated 
in a wounded nucleo n model of nuclei interactions in 
air (|Bialas et al.lll976D . 

5. COMPARISON OF FRAGMENTATION AND 
SUPERPOSITION MODELS 

In the current work, a simple eikonal approximation 
is used, implementing the Glauber model of nucleus- 
nucleus collisions. The geometric quantities — average 
impact parameter, number of wounded/spectator nucle- 
ons, etc. — are calculated using an algorithm de scribed 
by the PHOBOS collaboration (jAlver et al.ll2008t) . Only 
projectile nucleus fragmentation is considered, because 
the fragments of the target nuclei make an insignificant 
contribution to the cascade in the ultra-high energy do- 
main. In Fig. [5] the distributions of wounded-nucleon 
number as a ratio to the mass number are given for 
Fe and N projectile nuclei. Validity of the approxi- 
mation used is based on the correct description of D- 
Au collisions at 200 GeV in Relativistic Heavy Ion Col- 
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lider (jBialas et al.ll2005D . 

Inelastic cross sections for projectile nuclei with 
A 6 (1,56) wer e calculated by the Glauber method 
(jHeck et al.lll998f ). The next step is to calculate the ag- 
gregate shower by adding the sub-cascades induced by 
the primary nucleus fragments at different levels in the 
atmosphere. A simplification here is that the spectator 
nucleons of the projectile are assumed to subsequently in- 
teract as a nucleus with the mass number A — n w , while 
the wounded nucleons interact independently. For exam- 
ple, in the case of the primary nitrogen nucleus, we have 
a chain of N-B-Li-He-D-H nuclei interactions in air. 

In order to derive Ei and S^oo in EAS induced by 
fragmenting nuclei, we have to use the values as a func- 
tion of the initial nucleon energy and sub-shower initial 
points, Xi. In the case of the ionization integral, calcula- 
tions fo r three essent ially different models have demon- 
strated ()Ivanovll2007| ) a small variation of Ei (8Ei /Eq < 
0.007) in the energy interval (0.1 < E < 10) EeV and 
Xi € (0,200) g cm" 2 . 

Comparing the ionization integral in the fragmentation 
model 

k k 

and superposition model 

E? =AE t (^-,xl), 

we find a difference (Ef — Ef)/Eg in the two models of 
less than 0.5% for the mass number of primary nuclei in 
the interval (1,56). 

In contrast, there is a considerable difference in the 
charged-particle densities of the two models. While 5*600 
is almost proportional to the energy of primary pro- 
ton/iron nuclei according to the CO RSIK A/SIB YLL and 
CORSIKA/QGSJET simulations of lNagano et~all pOOOl ) 
(dlgE /dlgS^ Q ° = 1.015 ± 0.015), the density rises with 
Xj, at least for the arrays at sea level. iCapdevielle et al.l 
(20080 showed that A S6 oo ~ 400 g cm" 2 for models of 
high multiplicity, and As600 ~ 300 g cm -2 for models of 
low multiplicity with the primary energy Eq = 1 EeV. 

Calculation of Sgoo m the fragmentation and super- 
position models results in a 5g 00 /5| 00 ratio within the 
(1,1.17) interval for models of high multiplicity, and in 
(1,1.21) for models of low multiplicity. The relative dif- 
ference 5g 00 /5| 00 — 1 in the latter case, and [e[ —E-)/E , 
are shown in Fig. [3] as a function of the mass number of 
an EAS primary nucleus. 

Using a superposition model to derive the energy of the 
primary nucleus from the measured shower parameters 
Ei or 5*600, we find that E is correct with accuracy better 
than 0.5% for Ei measurement, and is overestimated by 
the factor Rs G (1,1-2) in the case of 5eoo- The actual 
values depend on the mass composition of UHECR flux. 

The measurements of the mass composition have been 
made by the HiRes collaboration in several energy bins 
above 0.1 E eV: i) HiRes prototype and the MIA muon 
array data (jAbu-Zavvad et abl 120011 ) led to the conclu- 
sion that CR intensity is changing from heavier to a 
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Fig. 4. — The differential energy spectra measured with EAS 
arrays. Energy corrections are applied assuming EAS induced by 
a primary iron nucleus. CR intensities and data symbols are the 
same as in Fig. [l] 

lighter composition between 0.1 EeV and 1 EeV; ii) then, 
Xmax distribution width and elong ation rate measure- 
ments in the interval (1,2.5) EeV (jAbbasi et al.l |2005|) 
are found to be consistent with a constant or slowly 
changing and predominantly light composition; iii) and 
finall y, the recent measur ements at energies above 1.6 
EeV (Abbas i et al. 2009b|) have found a proton domi- 
nated UHECR flux, within QGSJET01 and QGSJET-II 
two-component proton-iron model assumptions. 

On the contrary, PAO collaboration measurements of 
both the mean depth of shower maximum vs energy and 
RMS for the X max distribution suggest the co mposition 
becoming heavier above 3 EeV (|Cronin| [2009). In addi- 
tion, the muon data recorded by the Yakutsk array indi- 
cate the presence of a significant fraction of heavy nuclei 
in th e CR flux at energies E > 10 EeV (Glushko v et al.l 
2008): The proton fraction is estimated as 0.32 < fp < 
0.71 (95%CL) and iron fraction is 0.29 < f Fe < 0.68 
(95%CL) in the two-component (P+Fe) composition and 
EPOS interaction model. 

In this context we have to bear in mind both the pos- 
sibilities of light and heavy UHECR compositions. To 
illustrate the convergence/divergence of the energy spec- 
tra measured with EAS arrays, the surface array energies 
are decreased^ by the factor 1.2 (Fig. [4|) as a lower limit 
to the primary energy in the case of the iron nucleus as- 
sumed as a primary; other case of the proton primary is 
illustrated in Fig. [TJ In other words, two figures (Figs. [T] 
and [4| represent the two extreme cases of UHECR nu- 
clear compositions. 

6. RESULTS AND DISCUSSION 

It has been noted previously that UHECR en- 
ergy spectra measured with giant EAS arrays agree 
with each other if the energy scales a r e adjusted 
([Bahcall and Waxmanl (|2003fl ; rivanovl ()2007f ): lBerezinskvl 
(2008), and other papers cited therein). However, the 
energy correction factors needed to merge the spectra 
exceed experimental errors. For example, Bcrczinskv 
(2008) used values 1.2, 0.75, 0.83, and 0.625 to shift PAO, 



1 we have used here the attenuation length in sec 8 as an esti- 
mator of dxi/dlnSfjOO 



except HP data re-evaluated using COR- 
SIKA/QGSJET ijAve et al. 2003) 
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TABLE 2 

Energy correction factors for EAS arrays. 



log, (E ,eV) 

Fig. 5. — The observed spectra with energy scales adjusted. Cor- 
rection factors are taken from Table [2] CR intensities and data 
symbols are th e same as in Fig. Q] and [4] Model calc ulation re- 
sults: full line (Aloisio et al. 2007); d ot and dash lines IBerezhkd 
2009); dash-dot line (Wibig and Wolfcndalc 2005). 



AGASA, Akeno, and Yakutsk data to those of HiRes. 
Moreover, instrumental and modeling errors are believed 
to equiprobably increase or de crease the estimated en- 
ergy. Instead, as lWatsonl (|2009( ) concluded when compar- 
ing integral calibrated fluxes of UHECRs from giant ar- 
rays, Sqqq measurements and fluorescence measurements 
of EAS assemble in two separate groups of data (Volcano 
Ranch, Haverah Park, AGASA vs PAO, HiRes results). 

Considerations in this paper insist that, besides the 
'Moscow correction factorU which reduces all CR in- 
tensities measured with EAS arrays, there is another 
factor that reduces the estimated energy of primary nu- 
clei diminishing a difference between two groups of data 
(Fig. |4j . The remaining spread of E can be attributed 
to the real accuracy of the EAS measurement technique: 
instrumental errors + model uncertainty. 

Two variants of the energy correction factors for EAS 
arrays are given in Table [2] for the two extreme cases 
of the primary nuclei (H,Fe). Resultant spectra with 
adjusted energy scales are shown in Fig. [5] Energy cor- 
rection factors are assumed constant at energies above 1 
EeV, and the spectra are shifted to a sample mean. 

Having an indefinite mass of EAS primary nuclei and 
model uncertainties to assign the energy correction fac- 
tors, we can find an interval of values where Re would be. 
A sample of six experiments (Figs. HJ SJ estimates that 
the average energy determination error, 5E/E, inherent 
to giant EAS arrays, is within an (8,19)% interval. 

Now we can compare the measured energy spectra with 
model simulations of extragalactic (EG) CRs predomi- 
nating over the Galactic (G) component. The dip (de- 
scribed as 'bu mp' in early measurements) observed in the 
CR spectrum ([Dvakonov et al. 1978; Bo wer et al .11198 lh 
was explained by iBerezinskv a nd Grig or'eval (J1988|ll, 
who studied the Bethe-Heitler pair production of pro- 
tons from distant sources on the cosmic microwave back- 
ground. 



Array 


AGASA 


HP 


Yakutsk 


PAO 


HiRcs/TA 


r>b e 
E 


0.96 


0.98 


0.88 


1.15 


1.01 


nil 


0.81 


1.01 


0.71 


1.29 


1.09 



3 d erived by IZatsepinl [IT959T ). |M urzinl (|1965), and IKalmvkovl 
(fT969T > 

4 it has been update d a number of times since; the latest is made 
bv I Aloisio et al.l p007ft 



The authors conclude that the dip shape is formed 
by the universal modification factor and is independent 
of UHECR propagation details. The only phenomenon 
known to modify the dip is the presence of a significant 
fraction of nuclei in the primary beam. 

In Fig. [5] the spectru m calculated in the dip sce- 
nario (|Aloisio et al.ll2007f ) with protons accelerated in EG 
source^ is shown. Only the source luminosity is fitted 
to the experimental data. The position on the energy 
scale and shape of the dip agree with observed spectra 
within experimental errors, as well as the GZK effect. 
The excess-over- GZK flux observed by the AGASA ar- 
ray is considered to be the result of the primary-energy 
overestimation in incline d showers of the highest energies 
dCapdevielle et al.ll2008h 

IWibig and Wolfendald (2005) have argued that the an- 
kle in the spectrum marks the transition from G to EG 
components of CRs. This phenomenological 'ankle sce- 
nario' uses a sum of power law G and EG spectra with 
slopes differing by A7 ~ 1.8. The result is shown in 
Fig. [5] by the dash-dot line. The authors emphasize the 
sharpness of the ankle observed, d 2 In JE 3 /dln 2 E, which 
is consistent with the ankle scenario, while in the dip 
scenario the sharpness is of insufficient magnitude, espe- 
cially in the case of a large nuclei-to-proton ratio in the 
primary beam. 

The energy spectrum of UHECRs produced at the 
shock created by the expanding cocoons around ac- 
tive galactic nuclei combined with the G component of 
CRs produced in supernova remnants was calculated by 
iBerezhkol (|2009| ). Expected CR composition shows an 
increase of A at Eq ~ 0.1 EeV owing to the G compo- 
nent, and a second one at energy ~ 10 EeV, produced 
by nonrelativistic cocoon shocks. 

Calculated UHECR intensity (jBerezhkol [2009) as a 
function of energy is shown in Fig. [5] for the two cases: 
dip (dots) and ankle (dashed line) scenarios of the G to 
EG transition. Experimental errors are too large to dis- 
tinguish between alternative scenarios. 

7. SUMMARY 

A comparison is drawn between the energy spectra of 
UHECRs measured with giant EAS arrays, taking into 
account the necessary corrections to the measured CR 
intensity and the estimated energy of the primary par- 
ticle. The intensity correction factor for the measured 
CR flux due to a rapidly falling energy spectrum and 
instrumental errors was derived in Moscow many years 
ago. EAS modeling uncertainty in the primary energy 
estimation is due to unknown nucleus-nucleus interac- 
tion characteristics at energies far beyond those studied 
in accelerator experiments. 

5 acceleration spectrum index 7 = 2.7 
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In particular, there may be a difference between en- 
ergies estimated in the superposition and fragmentation 
models of nucleus-nucleus interactions. It is shown that, 
indeed, using the superposition model to analyze the sur- 
face array data, e.g. Seoo, we obtain E, which has to be 
corrected due to the fragmentation rate overestimated. 

Applying the essential correction factors, it is shown 
that UHECR energy spectra measured with EAS arrays 
are congruous within experimental errors arising from in- 
strumental and model uncertainties. Residual differences 
in the energy scales of giant arrays are used to estimate 
UHECR energy determination error inherent in EAS de- 
tection techniques. 
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